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^^ ■ Abstract. We implemented sink particles in the Adaptive Mesh Refinement (AMR) code 

Cn ' FLASH to model the gravitational collapse and accretion in turbulent molecular clouds and 

O , . cores. Sink particles are frequently used to measure properties of star formation in numerical 

^ ■ simulations, such as the star formation rate and efficiency, and the mass distribution of stars. 

^^ I We show that only using a density threshold for sink particle creation is insufficient in case of 

supersonic flows, because the density can exceed the threshold in strong shocks that do not 

necessarily lead to local collapse. Additional physical collapse indicators have to be considered. 

We apply our AMR sink particle module to the formation of a star cluster, and compare it 

to a Smoothed Particle Hydrodynamics (SPH) code with sink particles. Our comparison shows 

encouraging agreement of gas and sink particle properties between the AMR and SPH code. 
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^ ■ 1. Introduction 

^ ^^ ■ Stars form in t urbulent, magnetized mole cular clouds by local gravitational collapse of 

dense gas cores (JMac Low fc KlessenI 12004 ). Modeling this process in computer simula- 
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(S| ■ tions is extremely difficult. It is necessary to follow the freefall collapse of each individual 

^ \ star, while keeping track of the global evolution of the entire cloud at the same time. 

^ . The fundamental numerical difficulty is that the freefall timescale decreases with increas- 



ing gas density: tff = y37r732Gp. Modeling each individual collapse and following the 

Aq . large-scale evolution of the cloud over several global freefall times in a single magneto- 

hydrodynamical calculation is beyond the capabilities of modern numerical schemes and 

supercomputers. Thus, if one wants to model the evolution of such a cloud, the individual 

^^ I runaway collapse must be cut-off in a controlled way and replaced by a subgrid model. 

There are two such subgrid models: heating the gas or using sink particles. In the 

first approach the gas is heated up above a given densi ty threshold to preve nt artificial 

fragmentation beyond the resolution limit for collapse (jTruelove et al.l 119971 ). There are 



/\ I two problems with the heating approach. First, the Courant timestep decreases, because 
id . the sound speed increases, and second, the gas equation of state is changed above the 
density threshold. For molecular clouds, heating of the gas may only occur f or densities 



density tnresnold. l^or molecular clouds, neatmg oi tne gas may only occur t or densities 
p > 10~^^gcm~^ in the optically thick regime (|Larsonlll969l : |Penstonlll969 ). However 



gas can become denser than the threshold value in shocks that do not necessarily lead 
to the formation of a gravitationally bound structure. Thus, shocked gas not going into 
freefall collapse will be heated up artificially at least in the case where the threshold 
density for heating is in the isothermal regime. 

The alternative subgrid model is to use Lagrangian sink particles, a method invented 
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by lBate. Bonnell. fc Pried (|l995[ ) for Smoothed Partiele Hyd rodynamies (SPH), and first 
adopt ed for Eulerian, Adaptive Mesh Refinement (AMR) bv lKrumholz. McKee. fc KleinI 
(200J). If the gas reaches a given density, a Lagrangian, accreting sink particle is in- 



troduced. However, sink particles are supposed to represent bound objects, and thus, a 
density threshold for their creation is insufficient. Compression in shocks can temporarily 
create densities larger than the threshold without triggering gravitational collapse. Previ- 
ous grid-based implementations of sink particles are mostly based on a density threshold 
criterion. Her e, we present an im plementation of sink particles for the Eulerian, AMR 
code FLASH (JFrvxell et al.l 120001 ) that uses a series of checks, such that only bound and 



collapsing gas is turned into sink particles. We show that the star formation efficiency 
and the number of fragments is overestimated, if additional, physical checks in addition 
to the density threshold are ignored in the isothermal regime. 

The main results presented here were previously published in iFederrath et al.l (|2010h . 



However, we discuss two new tests for linear and angular momentum conservation in ^ 
In ^we present the physical checks necessary to avoid spurious sink particle creation 
and present the main results of the AMR-SPH comparison of sink particles. 



2. Sink particle creation checks and AMR-SPH comparison 

We refer the reader to IFederrath et al.l (|2010[ §2.2) for a detailed discussion and for 



the implementation of the sink particle creation checks. In summary, for successful sink 
creation the gas exceeding the density threshold must also 

(a) be converging (along each cardinal direction individually), 

(b) have a central gravitational potential minimum, 

(c) be Jeans-unstable (including magnetic pressure), 
{d) be bound (including magnetic energy), 

(e) be on the highest level of the AMR (i.e., Jeans length resolved), 

(f) not be within accretion radius of existing sinks (then accretion checks apply). 
During testing it turned out that the first two checks are particularly important to avoid 
spurious sink particle creation, however, the relative importance of each of the individual 
checks should be investigated in more detail (Wadsley et al., in prep.). 

Fig. [1] shows the comparison of FLASH (all checks on; top panel) against FLASH 
(first four checks switched off; bottom panel), clearly demonstrating the importance of 
the physical sink creation checks in the isothermal regime. The middle panel of Fig. [T] 
shows the SPH-NG run, exhibiting some differences to the FLASH run, which can be 
attributed to the slightly faster collapse of the cloud core in the SPH run. This is most 
likely a consequence of the faste r decay of the supportin g initial turbulence due to the 
slightly higher viscosity in SPH (jPrice fc Federrathll2010[ ). After correcting for this and 
comparing at times when 26% of the gas has been accreted onto sinks, the FLASH 
and SPH-NG runs are in very good agreement with 49 an d 50 sink particles for med, 
respectively, and having similar mass distributions (see §4 in IFederrath et al.ll2010l ). 

3. Momentum and angular momentum conservation test 

In IFederrath et al.l (|2010[ ) we performed a suite of test simulations for the new sink 



particle module in FLASH, including circular and highly eccentric orbits, the collapse 
of a Bonnor-Ebert sphere and a singular isothermal sphere, and a rotating cloud core 
fragmentation test. Here we add a momentum conservation test shown in Fig. |2l where 
we initialized a 0.46 Mq core (with very low resolution for testing purposes) at rest and a 
0.1 M0 sink particle with an initial momentum of Py = lO^^gcms"^ in ^/-direction. Both 
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Figure 1. The formation of a stellar cluster from 
freefall time. Top: FLAS H (default settings), m iddle. 



checks switched off). See iFederrath et al.l ()2010l ) for further details, 



100 M© turbulent cloud after one global 
SPH-NG, and bottom: FLASH (first four 



gas core and sink particle are offset from the center. Fig. [2] (left panels) show column 
density snapshots of the time evolution of the system. The initial gas core collapses and 
a second sink particle forms and accretes almost all the initial gas mass, while the initial 
sink particle only has time to accrete about 2% of the gas. The two sink particles are 
then followed for 20 orbits up to tend = 10^ yr. The right panels of Fig. [2] show the 
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log,„ ( column density [g/cm^] ) 




Figure 2. Left panels: time evolution of the column density in the momentum conservation 
test. Right panels: time evolution of the momenta of sink particles and gas in x- (top) and 
y-direction (bottom). For momentum conservation to hold, the momenta of sinks and gas must 
be symmetric about the dotted line, which indicates half the initial total momentum. 
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Figure 3. Left: Column density snapshot of a variant of the Boss fc Bodenh eimer (1979") test, 
followed to very late times when 12 fragments have formed and 68.5% of the gas has been 
accreted by sinks. Right: Shows that angular momentum is conserved to within 0.5% for all 
times. 



time evolution of the momentum. Momentum is conserved to within 3% over the whole 
duration even in this very low resolution run, shown by the symmetry of sink and gas 
momenta. Kinks in the momentum evolution indicate strong accretion events. 

For testing angular momentum conservation, we use a high-resolution model of the 
IBoss fc Bodenheimer! (|l979[ ) test. Figure [3] (left) shows the column density of the disk 
(face on in the x-?/-plane) for the last time frame computed, where 12 fragments have 
formed in the disk, which turned into a ring-like structure at late times. The right panel 
shows the time evolution of the total angular momentum, L^. Angular momentum is 
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Figure 4. Left: Density slice through a rotating, magnetized disk. A bipolar outflow (velocity 
vectors are overlaid on the density) has formed and propagated roughly 5, 000 AU from the disk 
midplane. The central sink particle has a radius of 80 AU. Right: Shows the three-dimensional 
magnetic field structure in the inner 800 AU of the disk. The magnetic field is wound up strongly 
and launches the outflows shown in the left image. 



conserved to within 0.5% (see ICommercon et al.ll2008L for a comparison of angular mo- 
mentum conservation in SPH and AMR). 



4. Conclusions 

We introduced and tes ted a new sink particle method for the AMR code FLASH 



Federrath et al.l (|2010[ ). In addition to the tests shown there we presented a linear 



and angular momentum conservation test in ^ A comparison of g;a s and sink particl e 
properties showed encouraging agreement with the SPH-NG code (JBate et al.l 119951 ) . 



More recently, another SPH code, GASOLINE (IWadslev et al.ll2004I ). also showed very 
good agreement, if the sink particle creation checks outlined in ^are used (Wadsley et 
al., in prep.). 

Some open issues remain concerning the modeling of ma gnetic fields in combination 
with sink particles in SPH (|Price fc Bate! 120081 : |Pricell2010[ ). The problem is that SPH 
particles are accreted inside the sink particle radius and are thus lost as resolution el- 
ements for the magnetic field. The advantage of using a grid-based implementation of 
sink particles is that the stencil for the magnetic field remains when gas is accreted. The 
geometry of the magnetic field thus remains intact. We also performed simulations of 
the collapse of magnetized, rotating cloud cores, which self-consistently produce bipolar 
outflows (Duffin et al.; Seifried et al., in prep.), showing that our sink particle approach 
works in combination with magnetic fields (Figure H]). 
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